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Abstract 

We present a 'new generation' model for high energy proton-proton 'soft' interactions. 
It allows for a full set of multi-Pomeron vertices, as well as including multichannel eikonal 
scattering. It describes the behaviour of the proton-proton total, atot, and elastic dae\/dt, 
cross sections together with those for low and high-mass proton dissociation. Although 
the model contains a comprehensive set of multi-Pomeron diagrams, it has a simple par- 
tonic interpretation. Including the more complicated multi-Pomeron vertices reduces the 
absorptive effects as compared to the predictions in which only the triple-Pomeron ver- 
tex is considered. Tuning the model to describe the available 'soft' data in the CERN 
ISR - Tevatron energy range, we predict the total, elastic, single- and double-diffractive 
dissociation cross sections at the LHC energy. An inescapable consequence of including 
multichannel eikonal and multi-Pomeron effects is that the total cross section is expected 
to be lower than before: indeed, we find cJtot — 90 mb at the LHC energy. We also present 
differential forms of the cross sections. In addition we calculate soft diffr active central 
production. 



1 Motivation 



It is essential to have a good model for the soft interactions of hadrons at high energies to, 
in particular, predict at the LHC, (i) the structure of underlying events, (ii) the value of 



the total cross section, (Jtot and the behaviour of the elastic cross section, da^i/dt, and (iii) 
the probability of diffractive dissociation. Moreover, it is extremely important to understand 
the asymptotic behaviour of high energy interactions in a partonic framework. This will give 
the possibility of simulating the underlying events using a Monte Carlo generator based on a 
theoretically justified description of soft interactions. In addition, we need such a model to 
account for the effects of the underlying event in high energy hard interactions and to calculate 
the survival factors of rapidity gaps in exclusive and other hard diffractive processes. Finally 
it is crucial to have a good understanding of the spectra of leading nucleons, that is of the 
diffractive dissociation cross section, in order to describe extensive air showers and to interpret 
the highest energy cosmic rays. 

2 Overview of existing descriptions of soft interactions 

First we have the very simple Donnachie-Landshoff parametrization [1] of the high energy elastic 
amplitude, which is described in terms of simple poles in the complex angular momentum plane: 
namely the Pomeron and secondary Regge poles. The latter give a negligible contribution at 
Tevatron and LHC energies. However this parametrization says nothing about the distribution 
of secondary particles in the underlying event. Moreover, already at the LHC energy, ^/s — 14 
TeV, the amplitude A{b, s) violates the black disc unitarity limit at small impact parameters 

We need to satisfy, at least, the two-particle s-channel unitarity relation, in order to respect 
the Froissart bound, and to describe the elastic cross section, which is about 20-25% of the 
total cross section in the Tevatron to LHC energy range. This leads to the eikonal form of 
the elastic amplitude, see (6) in Section 3. To allow for the possibility of proton excitations 
we need to consider multichannel (say n-channel) eikonal models which include rescattering in 
the i = I, ...n diffractive eigenstates [2, 3, 4]. We will review the eikonal approach in the next 
Section. At present one- and two-channel eikonal models are used to predict the LHC cross 
sections [5, 6, 7]. A one-channel approach was used in [6, 8], and two-channel eikonals were 
used in [7, 5, 9]. 

However, even a multichannel eikonal is unable to account for diffractive dissociation into 
high-mass states. These processes are usually described in terms of Regge theory with the 
help of the triple- Pomeron vertex. The problem is that, after we allow for the low values of 
the probability, S^, that the rapidity gaps survive the eikonal rescattering, the value of the 
triple- Pomeron coupling, g^jp, needed to describe the data (a ~ S^g^p), becomes rather large; 
namely g^jp ~ 5'Ar/3 [10], where gi\f is the nucleon-Pomeron coupling^. As a result we cannot 
neglect more complicated multi- Pomeron diagrams containing a large number of triple- Pomeron 

^Earlier estimates [11, 12, 3], which do not account for the rescattering factor S*^, give g^p ~ gAr/10. 
However, as discussed in [13, 14], when we account for more complicated enhanced diagrams we need a larger 
value of the bare triple- Pomeron coupling. 
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Figure 1: A typical multi-Pomeron fan diagram contributing to the interaction between a small 
size beam particle and a large size target particle. 

couplings. A specific set of such diagrams, known as "fan" diagrams, have been summed in 
Refs. [15, 16, 17]. The selection of these fan diagrams is justified for proton-nuclei interactions 
or deep inelastic scattering. That is when the size of the incoming object is much smaller than 
the size of the "target" . In other words, we account for the multiple interactions with the large 
size "target", but not with the "beam" particle, see Fig. 1. The summation of these diagrams 
leads to saturation, that is the amplitude T constant at large s. The problem is that 
the fan diagram approximation can only be used to describe the beginning of the approach 
to saturation. Once we are near saturation it is not justified to neglect more complicated 
multi-Pomeron graphs^. 

A good description of CERN ISR-Tevatron data on atot, daei/dt was obtained in [7], where 
the high-mass diffractive dissociation was described phenomenologically and added to the two- 
channel eikonal. That is, high-mass diffraction was not generated from the underlying theory, 
but was included in terms of the leading triple-Pomeron diagram, and then this contribution 
was added to the total proton opacity Q{b). 

The shortcomings of the existing approaches to describe soft interactions can be summarised 
as follows. First, it is not clear, a priori, how the results will change in going from a two- to a 
three- (or more) channel eikonal model. Furthermore, we do not know how the results depend 
on the size (that is, the form factor) of each diffractive eigenstate^ i. Moreover, since g^jp is 
not small, we cannot neglect the more complicated multi-Pomeron interactions^. 

Here we will consider a three-channel eikonal model. We will assume that the cross section 
of each eigenstate, cxj (that is, the coupling of the eigenstate to the Pomeron) is proportional 
to the square of the transverse size of the component, R^. This form is motivated, either 
by leading-order QCD, where a ~ otsR^-, or by the assumption that each eigenstate has an 

more general, but still incomplete, set of multi-Pomeron graphs generated by g'^p was considered in 
Ref. [18], and more recently in Refs. [19, 20]; note that in [19], besides (?3ip, more complicated multi-Pomeron 
vertices were considered 

^In [7] both eigenstates are assumed to have the same size, and in the latest versions of the Tel- Aviv model 
[9] another extreme is considered — the size of the second component is zero. 

^Here we seek a model for the high energy pp interaction and therefore we neglect secondary Reggeon 
contributions, and include only the Pomeron. 
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opacity of the same shape with the same value of Qi{b — 0), such that the integrated couphng 
/ (Pb Jlj oc i?^. Another possibihty, which we will analyse, is to use the same form factor, that 
is the same size, for each component, as in [7]. 

By allowing for a different size for each component, we have the possibility to introduce, 
and to calculate the different partonic composition of each diffractive eigenstate. For instance, 
for a smaller size component, the evolution will start at a larger scale and so there will be 
less gluons and other low x partons radiated during the evolution. 



Let us delay the inclusion of high-mass diffractive dissociation for the moment in order to 
first introduce the eikonal description of soft proton-proton interactions. This will allow low- 
mass diffractive excitations to be included. We will then have to describe how to incorporate 
the important contributions made by high-mass single and double proton dissociation. The 
high-mass dissociations have a crucial effect on the predictions for the LHC. 

Unitarity plays a pivotal role in diffractive processes. The total cross section is intimately 
related to the elastic scattering amplitude and the scattering into inelastic final states via 
s-channel unitarity, SS^ = /, or 



with S — I + iT. 

3.1 Single- channel eikonal model 

First, we briefly recall the relevant features of the single-channel eikonal model. That is we 
focus on elastic unitarity. Then "disc T" is simply the discontinuity of T across the two-particle 
s-channel cut. At high energies we have a sizeable inelastic component. The s-channel unitarity 
relation is diagonal in the impact parameter, 6, basis, and may be written 



3 Eikonal model for diffractive scattering 



disc T = T-T^ ^ iT^T 



(1) 



2Imrel(s, b) = |re,(s, 6)1' + Gine^S, b) 



(2) 



with 




(5) 



(4) 



(3) 
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These equations are satisfied by 



ImTei(s,6) = (6) 
aei(.,6) = (l-c-^/^)^ (7) 
ainei(s,6) = l-e"^, (8) 

where Q{s,b) > is called the opacity (optical density) or eikonal^. Prom (8), we see that 
exp(— Q(s, b)) is the probability that no inelastic scattering occurs. 

The well known example of scattering by a black disc, with ImTei — 1 ior b < R, gives 
Cei — Cinei = Tfi?^ and (Ttot = 2t:B?. In general, we see that the absorption of the initial wave 
due to the existence of many inelastic channels leads, via s-channel unitarity, to diffractive 
dissociation. 



3.2 Inclusion of low- mass diffractive dissociation 

So much for elastic diffraction. Now we turn to inelastic diffraction, which is a consequence 
of the internal structure of hadrons. This is simplest to describe at high energies, where the 
lifetime of the fluctuations of a fast hadron is large, r ~ £'/m^, and during these time intervals 
the corresponding Fock states can be considered as 'frozen'. Each hadronic constituent can 
undergo scattering and thus destroy the coherence of the fluctuations. As a consequence, the 
outgoing superposition of states will be different from the incident particle, and will most likely 
contain multiparticle states, so we will have inelastic, as well as elastic, diffraction. 

To discuss inelastic diffraction, it is convenient to follow Good and Walker [2], and to 
introduce states 0^ which diagonalize the T matrix. Such eigenstates only undergo elastic 
scattering. Since there are no off-diagonal transitions, 

(0,|r|0,) = O iorj^K (9) 

a state k cannot diffractively dissociate into a state j. Wc have noted that this is not, in general, 
true for hadronic states, which are not eigenstates of the ^'-matrix, that is of T. To account for 
the internal structure of the hadronic states, we have to enlarge the set of intermediate states, 
from just the single elastic channel, and to introduce a multichannel eikonal. We will consider 
such an example below, but flrst let us express the cross section in terms of the probability 
amplitudes of the hadronic process proceeding via the various diffractive eigenstates^ ^jfc. 

Let us denote the orthogonal matrix which diagonalizes Im T by a, so that 

ImT = aFa^ with {(pj\F\(pk) = Ft 6jk. (10) 

■''Soinctimcs 0/2 is called the eikonal; for simplicity we omit below the real part of Tei. At high energies, the 
ratio ReTei/ImTei is small, and can be evaluated via a dispersion relation. 

^The exponent exp(— fifc) describes the probability that the diffractive eigenstate (pk is not absorbed in the 
interaction. Later we will see that the rapidity gap survival factors, S"^, can be described in terms of such 
eikonal exponents. 
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Now consider the diffractive dissociation of an arbitrary incoming state 



k 

The elastic scattering amphtude for this state satisfies 

(illmTli) = Yl M'Fk = (F), (12) 

k 

where = and where the brackets of (F) mean that we take the average of F over 

the initial probability distribution of diffractive eigenstates. After the diffractive scattering 
described by Tfj, the final state |/) will, in general, be a different superposition of eigenstates 
from that of which was shown in (11). At high energies we may neglect the real parts of 
the diffractive amplitudes. Then, for cross sections at a given impact parameter b, we have 



d<j. 



tot 



d% 

d(Te\ + ST> 



= 2Im(j|T|j) = 2 ^ \a,k?Fk = 2(F) 

k 

- mm? = (e mf^' = {F? (13) 



d% 

It follows that the cross section for the single diffractive dissociation of a proton, 

^ = {F') - {Ff, (14) 

is given by the statistical dispersion in the absorption probabilities of the diffractive eigenstates. 
Here the average is taken over the components k of the incoming proton which dissociates. If 
the averages are taken over the components of both of the incoming particles, then in (14) we 
must introduce a second index on F, that is Fn^, and sum over k and i. In this case the sum is 
the cross section for single and double dissociation. 

Note that if all the components (f)k of the incoming diffractive state |j) were absorbed equally 
then the diffracted superposition would be proportional to the incident one and the inelastic 
diffraction would be zero. Thus if, at very high energies, the amplitudes Fk at small impact 
parameters are equal to the black disk limit, F^ = I, then diffractive production will be equal 
to zero in this impact parameter domain and so will only occur in the peripheral b region. A 
similar behaviour already takes place in pp (and pp) interactions at Tevatron energies. Hence 
the impact parameter structure of inelastic and elastic diffraction is drastically different in the 
presence of strong s-channel unitarity effects. The elastic amplitude originates mainly from the 
centre of the disk (that is, from small b), while dissociation comes from the periphery. Hence 
it is important to pay special attention to the periphery of the proton, in impact parameter, 
b, space. First, large values of b are responsible for the small t behaviour of the amplitude. 
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Second, the large b region, where the optical density (or opacity), becomes small, gives 

the major contribution to the survival probability of rapidity gaps. 

It is clear that if we allow for the possibility of diffractive excitation, then we will enlarge 
the absorptive effect caused by the eikonal. As a consequence the results obtained by a one- 
and two-channel eikonal fit will be substantially different. However, once we fix the dispersion, 
(14), that is, once we fix the ratio''' o"sd (low- mass) /cd, it turns out that the inclusion of a third 
channel will not practically change the result. Of course, with a three-channel eikonal we have 
too many parameters. We will therefore consider an extreme case where each channel has the 
same weight; = l/\/3 with i = 1,2,3.. Otherwise the contribution of the eigenstate with 
the smallest weight will be less visible and the situation will be close to the 2-channel model. 
Next, we fix the value of the dispersion, (14), 



As a result the total probability of dissociation (in the limit of small opacity, Q <^ 1), 



will be the same in the two- and three-channel cases. Here (3i is the 0j-Pomeron coupling, so 
the amplitude Fj is proportional to the factor of 2 in (16) accounts for low-mass excitations 
of both the colliding beam particles. 

Note that to describe the elastic scattering data in terms of a multichannel eikonal model, 
we have to increase the size of the intercept, q;jp(0), and to reduce the slope a', of the bare 
Pomeron trajectory, in comparison with the naive Donnachie-Landshoff amplitude, in order to 
compensate for the absorptive effects of multichannel eikonal rescattering [3, 14, 21]. 

It was checked that the effects of the pion loop in the bare Pomeron trajectory (see [7]) 
can be mimicked by increasing the value of a'. Thus, if we do not want 1% accuracy, we 
can replace our previous two-channel description [7] of the data by a simpler two-channel fit, 
which is suitable for the calculation of the rapidity gap survival factors for the different hard 
diffractive processes and to evaluate the soft cross sections atot, daei/dt and crgg"^. 

3.3 Survival factors 

An essential ingredient, in the prediction of the rate of a particular diffractive process, is the 
calculation of the suppression due to multi-Pomeron exchanges^. We assume that the rapidity 
gap is large enough that it will be populated by secondaries from soft rescattering; that is, in an 
inelastic soft interaction, there is a negligible probability for a fluctuation with zero multiplicity 

^Note that in terms of a limited number of eigenstates (pi, we can consider only 'low-mass' proton excitations, 
within a limited interval of M^. 

^This is especially important for searches for New Physics signals in processes with tagged forward protons, 
see, for example, [22, 23]. 




(15) 



_lowM / _ Q. 2 



(16) 
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within the large gap interval. Then we can write the survival factor of the gap with respect 
to the i — k soft interaction at fixed impact parameter b as e'xp(—Qik{b)). This survival factor 
should be averaged over the full set of diffractive eigenstates and over the impact parameter b. 
If we do not detect the outgoing protons the survival factor caused by eikonal rescattering is 
given by 

jd%\api\'^ |ap/fep lA^ifeP exp(-l]ifc(s, 6)) 
= ^ J , (17) 

i,k 

where |apip and |ap'fcp are the probabilities of finding the partonic diffractive eigenstates i = 
and k = 10^) in the two colliding proton states \p) and \p'), respectively, see (11). |A^tfcP is the 
probability of producing the particular final system from the incoming eigenstates i and k. If 
one, or both, of the outgoing protons is observed, then we have to average the amplitude and 
not the cross section. For example when both protons are tagged with = we have 




The survival probability, S'^, can depend on the types of active partons j and the value 

— * 

of their momenta, Xj and k±j. In a multichannel eikonal model, with a few eigenchannels, 
we expect the channel with the smallest cross section to contain mainly valence quarks with 
larger x, while the channels with larger cross sections will be due to sea quarks and gluons, 
concentrated at smaller values of x [24] . 



4 Inclusion of high-mass diffraction 

In the eikonal approximation the interaction of two diffractive eigenstates, (pi and (pk was 
described by the opacity fli^ corresponding to the two-particle s-channel irreducible amplitude. 
This flik generates the complete two-particle reducible amplitude Fik, 

F,k{s,b) = l-e-^-/^ (19) 

see (6). The connection between the irreducible amplitude fik and the full amplitude Fik is 
illustrated in Fig. 2. The defect of the approach, where the irreducible amplitude Qik = fik of 
the H — k^ interaction has been approximated by a single pole in the complex angular momentum 
plane, is that we are unable to separate the contribution coming from events with a rapidity 
gap from those of completely inelastic multiparticle processes. Thus we are unable to identify 
high-mass diffractive dissociation. Moreover, if we try to describe high-mass dissociation in 
the Good- Walker formalism we encounter the problem of double counting when the partons 
originating from the dissociation of the beam and 'target' initial protons overlap in rapidities. 
These problems are well illustrated by the soft interaction sketched in Fig. 3. In this Section 
we discuss how to treat high-mass diffractive dissociation. 
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Figure 2: The connection between the full amplitude Fik and the irreducible amplitude fik for 
scattering between i and k diffractive eigenstates. 




Figure 3: A 'soft' high energy interaction in which partons originating from the dissociations of 
the colliding protons overlap in rapidity. The overlap illustrates the impossibility of describing 
high-mass diffraction in terms of a pure eikonal (Good- Walker) formalism. Note that the master 
equations which describe the evolution of these 'parton showers' are introduced in Section 4.2. 
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Figure 4: The addition of the triple- Pomeron diagram which allows for high-mass, M, diffractive 
dissociation. 
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Figure 5: The diagrams on the left-hand-side constitute the Schwimmer model. We also show 
a more complicated diagram. 



4.1 The deficiency of the Schwimmer model 

The usual way to include large mass M dissociation is to include the triple-Pomeron diagram 
as in Fig. 4. However, the triple-Pomeron couphng is not small, g^p ~ gf^r/S [10], so to be 
consistent we have to sum up a series of more complicated graphs like those in Fig. 5. As 
mentioned above, if the size of one object "fc" is much larger than that of the other, "i", 
then the major contribution comes from the sum of fan diagrams with multiple interactions on 
k. This is the so-called Schwimmer model [15]. In the Schwimmer model^ the cross section 
becomes saturated at high energy, that is at high Y, [15, 25, 26] 

9i9k e^^ _ giQk , . 

where Y = Ins is the available rapidity interval and A = ajp(O) — 1; and where the value of 
e = gsipgk/^ fixes the maximum saturation density that can be reached in the Schwimmer 
model. At the same time, large mass diffractive dissociation is given by 

, largeM a A(2Y~lnM^) , largeAf n n A 

?.^2 "^SD ^ 9i9k A e e ^ Af^ = ^^^^ foil 

dM^ [l + e(2eAi'_e^(^-i'-A'^') -1)]2 dM^ 4e(M2)^' ^ ^ 

where the final expression applies asymptotically, when 2e e^^ ^ 1. By comparing (20) and 
(21), we see that the absorptive effects are much stronger in diffractive dissociation than in 



^For simplicity, we use the model in 1+1 space-time dimensions. We therefore use the notation g^p for the 
triple-Pomeron coupling in 1+1 dimensions in order to distinguish it from the real triple-Pomeron coupling gsp. 
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the total cross section. Note the second power in the denominator in (21) and, moreover the 
denominator contains 2exp(Ay) rather than exp(Ay). This reflects the well known fact that 
the screening of the diffractive dissociation amplitude is twice as strong as that for the elastic 
amplitude. Recall that eqs. (20) and (21), and the discussion in the following subsections, 4.2 
and 4.4, are in 1+1 dimensions. Non-zero transverse momentum is introduced in Section 5. 

If, in the real world, we try to describe the total cross section in terms of the Schwimmer 
model we obtain a diffractive cross section, da^D / dM"^ , which is an order of magnitude lower 
than that observed at the Tevatron. Indeed, as was discussed in Section 2, the triple-Pomeron 
vertex cannot be too small. Even neglecting the screening corrections, we need g^jp ^0.1. 
In particular, with g^ip/gN = 0.1 — 0.2 we flnd-'^^ that the cross section dasu/dxL is about 5 
times smaller than that observed by the CDF collaboration at the Tevatron. With a larger 
gsp the "e" term in the denominator of (21) already becomes large for Tevatron energies; 0(1) 
for g^p > 0.5. Moreover, from the rightmost 'limiting' expression in (21), we see that aso is 
smaller for larger values of g^jp. For larger g^jp = gN/'2 we found that saturation is quickly 
obtained. However, as seen in Fig. 6, it is still possible to describe the total and elastic cross 
section data. Here, to compensate for the strong absorption, we need A = 0.5. At the Tevatron 
energy^^ this results in a diffractive cross section 

(1-^l)4^ = 0.2mb/GeV^ for i = 0, = 0.99, (22) 
dxidt 

which is orders of magnitude less than that measured by CDF [27, 28]. 

Next, dasT)/dxL decreases with energy as s""^, since for a fixed longitudinal beam mo- 
mentum fraction, xl, carried by the recoil proton^^, we have = (1 — xl)s. Moreover, as 
discussed in [29], the strong absorptive effect implied by the fan diagrams is not observed in 
the leading neutron spectra obtained at HERA. So the Schwimmer model is in contradiction 
with experiment. 



4.2 Including more complicated multi-Pomeron diagrams— zero trans- 
verse dimension 

For simplicity, we continue the discussion in 1+1 dimensions and wait until Section 5 to intro- 
duce non-zero transverse momentum. It is important to note that, besides the triple-Pomeron 
vertex, there are more complicated multi-Pomeron interactions, see Fig. 7. Indeed, it was 
proposed [30, 31] to sum diagrams with all possible (n, m) configurations at each vertex with 
coupling gl^. The exact values of the various vertices are, of course, unknown. In Ref. [30] 

^°The procedure that we use to include the transverse size of the proton, and other practical details, are 
described in Sections 5 and 6. 

^^The Tevatron data that we use here were collected at an energy ^/s = 1.8 TeV. For this reason we use a 
Tevatron energy of 1.8 TeV when we present our predictions. 

^^That is, for a fixed rapidity gap, Ar] = ln(l/(l — xl))- 
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Figure 7: A multi-Pomeron vertex 



11 



Figure 8: A typical parton shower. 



Cardy proposed a beautiful way to study the asymptotic high energy limit, by assuming the 
analyticity of (7^ in the Re m, Re n > semi-plane. However, to describe the cross section 
at the energies available to experiment we need to know the explicit form of the vertices. The 
simplest possibility is to assume an eikonal form 

gl cx X^^^-\ (23) 

Such a programme was carried out in 1986 [14], and led to a rather reasonable prediction at 
the LHC energy. 

Here we will consider the partonic interpretation of an analogous approach. As mentioned 
above, for simplicity, we first neglect the transverse dimensions. In terms of the partonic 
evolution in rapidity space, the bare Pomeron pole contribution to the elastic amplitude / is 
generated by the simple equation 

= A /(,). (24) 

This equation can be regarded either as the equation for the single Pomeron amplitude, or 
as the equation for the parton density generated from either the beam i or the target k. For 
definiteness, let us consider the partonic evolution from the target k. Then (24) gives 

f = 9l (25) 

where A corresponds to the intercept of the pole amplitude. That is^^, / oc with a = 1+A, 
so that a = gl{s/so)^, where gk is fixed by the initial condition, namely the probability of the 
interaction at s = Sq — in other words, by the parton density at y = 0, the beginning of the 
evolution. 

Here, at a low scale, it is impossible to say whether a parton is quark or gluon. These 
degrees of freedom are not well defined in the 'soft' regime^^. However, in a soft high energy 
interaction, the parton may be considered as a 'small' elementary object which mediates the 
process, in the spirit of the original parton model [32, 33]. Thus, here, by the partonic picture 

^^Here, for fixed impact parameter, we use the same normalization of the amplitude as in (3). 
^^Nevertheless we hope that at larger scales our partons will match smoothly with quark and gluons. 



12 




Figure 9: The evolution of the elastic bare Pomeron amplitude 
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Figure 10: The ladder structure of the triple-Pomeron amplitude. 



we mean the evolution of the parton shower as shown in Fig. 8. We emphasize, that our 
approach goes beyond simply the probabilistic partonic interpretation of the amplitude, but, 
more important, it allows us to describe the actual evolution of the parton shower through 
master equations of the form of (24) (or (26), (29) and (30) below). 

Indeed, according to (24), the 'splitting function' P, which specifies the probability of the 
emission of an extra parton c within the rapidity interval dy, is simply P = A. When we iterate 
(24) we generate a ladder-type amplitude, which has exactly the structure of that considered 
in Ref. [34], but now in 1+1 dimensions. This is illustrated in Fig. 9, which produces an 
exchanged Pomeron (shown by the bold dashed line) with a ladder-type structure. 

A multi-Pomeron contribution arises from the absorption of the intermediate s-channel 
partons c during the evolution of / in y. In particular the triple-Pomeron diagram in Fig. 4 
means that parton c undergoes an extra rescattering with the target parton k, as shown in 
Fig. 10, [33]. Allowing for many rescatterings, we have to sum over different numbers of ladders 
between partons c and k. Assuming an eikonal form for the multi-Pomeron-proton vertex, it is 
natural to replace (24) by 

f = fA exp(-fi/2), (26) 
dy 

where the 'opacity' Q = Xf describes the transparency^^ of the target k. Since we are dealing 

^^Since the opacity depends on the type of mcoming particle, it is natural to introduce a parameter A which 
reflects the difference of the opacity of the target felt by the intermediate parton and the opacity felt by the 
incoming eigenstate i. 
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with the elastic amplitude f we use exp(— Q/2) and not exp(— Q), see eq. (3). The coefficient 
A reflects the fact that parton c may be different, and have a different absorptive cross section, 
from that of the cigenstate i. The value of this parameter A should be tuned to describe the 
data driven by the triple-Pomeron vertex. Since the triple-Pomeron vertex is relatively small, 
we expect A < 1. Indeed, as we have already discussed, the HERA data on inelastic J/ip 
production (with target proton dissociation) indicate that [10] Qsip/qn ~ 0.3. In the limit 
A — > 0, we come back to the pole amplitude of (24). For non-zero A, it is straightforward to 
check that the asymptotic amphtude generated by (26) grows slowly as 

f ^ J [In(yA) + In \ny + ....] . (27) 

This is very different to the Schwimmer model, which only accounts for the triple-Pomeron 
vertex. There we kept only the first two terms in the decomposition of the exponent in (26), 
so that the equation took the form 

f ^ Af - gsjp f. (28) 
dy 

In that case the value of / rapidly saturated, / — > constant; see (20) with (Jtot ~ /■ The crucial 
difference is that in the perturbative (with respect to small g^p) calculation, the right-hand- 
side of (28) may even be negative for large /; while, after resummation of the multi-Pomeron 
exchanges, the right-hand-side of (26) is definitely positive 

In terms of Regge diagrams, (26) sums up the system of fan diagrams in which any number 
m of "lower" Pomerons couples to a fan vertex g^^, defined as in Fig. 7. In order to include the 
rescattering with the beam i we replace (26) by 



dy 



= f{y) A e-(''^^y)+''^(y >^/\ (29) 



The final term in the exponent is the opacity of the beam i, which depends on the rapidity 
interval y' — Y — y, with Y — Ins. The equation for the opacity fli = Xfi has the analogous 
form ^ 

^ = f(y') A e-("'(2'')+"*(f))/2. (30) 



dy' 

Actually, this is the same equation as (29), but now evolving in the backward direction starting 
from the boundary condition f{y' = 0) = gi at y = Y . 

The system of equations (29) and (30), with boundary conditions f{y — 0) — gk and 
f{y' = 0) = gi may be solved by iteration [19, 20]. Depending on the values of gi, gk and A, we 
usually need no more than 5 — 15 iterations to reach an accuracy of 0.1%. That is the forward 
evolution of the amplitude f{y) in the 'background' field f{y') gives the same result, to 0.1% 
accuracy, as the backward evolution of f{y') in the 'background' field f{y)- Once (29) and (30) 
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For the case of a nuclear target k, an equation analogous to (26) was discussed in [35]. 



14 




Figure 11: Symbolic representation of the sum of multi-Pomeron diagrams generated by (a) 
equation (26), (b) equation (29), (c) equation (30), and (d) the system of equations (29) and 
(30). 



are solved, it is straightforward to simulate in a Monte Carlo the development of the parton 
shower (29) in the known external field Qi. 

At this stage, it is useful to have a pictorial illustration of the above equations. In terms 
of the diagrams of Reggeon Field Theory, (28) produces the set of fan diagrams generated by 
the triple-Pomeron vertex only. The first two diagrams are shown in the left part of Fig. 5. 
On the other hand, (26) produces an analogous set of fan diagrams, but now the number of 
"lower" Pomerons in any vertex is arbitrary. Expanding the exponent on the right-hand-side 
of (26) leads to vertices gl^ with m > 2. The factor 1/m! which comes from the expansion 
of the exponent accounts for the identity of the bare Pomerons^^. The set of fan diagrams 
generated by (26) is shown symbolically in Fig. 11(a) as a shaded triangle. Now we turn to 
equations (29) and (30). The factor exp(— fij(?/')/2) in (29), which describes the rescattering 
of the intermediate parton c with the beam i, generates vertices like g^, and leads to more 
complicated diagrams, one of which, for example, is shown by the last diagram in Fig. 5. If we 
expand both exponents exp(— nfc(j/)/2) and exp{—fli{y')/2), then we obtain diagrams with a 
full set of vertices gl^. In Fig. 11(b), the effect of the absorptive factor exp{—fli{y')/2) is shown 
by a dashed line. The opacity f2j is generated by (30), which corresponds to Fig. 11(c). Solving 
the system of equations (29) and (30) we obtain the full amplitude, which we show symbolically 
as the light-shaded rectangle in Fig. 11(d). 

Note that the fan diagram of Fig. 11(a) contains a contribution with a rapidity gap in 
some interval {0,y), as shown in Fig. 12(a). An analogous rapidity gap can be found in the 
whole amplitude Fig. 11(d). We show this contribution in Fig. 12(d), where the light-shaded 
'rectangles' correspond to total amplitudes generated by (29) and (30). Note that in each 
'rectangle' in Fig. 12(d) we account for the absorptive effect exp(— fi/2) of the corresponding 
external field generated in the whole rapidity interval {0,Y), and not just in the sub-intervals 
{0,y) or {y,Y) occupied by the particular 'rectangle'. 

^^Strictly speaking, equations (26), (29) and (30) correspond to vertices with g'^ cx nmA"+™~^. 
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Figure 12: The symbolic description of diagrams appropriate for diffractive dissociation of the 
eigenstate i with a rapidity gap in some interval {0,y). The notation (a) (b,c) (d) refers 
to an analogous development to that indicated in Fig. 11. The cross on the lower line indicates 
that eigenstate k is on-shell. 

4.3 Comments on the final state 

Recall that in the evolution equations for amplitude, given in (26), (29) and (30), we include the 
absorptive factor exp(— fi/2) and not exp(— f2). That is we work with the forward amplitude 
ImT = 1 — e~^/^, which at each step of the evolution (in rapidity y) includes all possible 
processes - both elastic and inelastic interactions with cross sections cxei = (1 — e~^^'^)'^ and 
CTinei = 1 - e"^; where cxtot = 2ImT = cXei + Uinei, see eqs. (6)-(8). 

Strictly speaking using the AGK cutting rules [36], together with vertices gl^ oc nmA"^™""^, 
we would obtain a parton-target elastic cross section cjei = (1 — e^^/^)(i7/2)e^^/^, instead of 
form (7). On the other hand the AGK rules were not proved for the general case of gl^ with 
n + m > 3. Here we prefer to adopt the partonic basis and to calculate the elastic cross section 
following the usual relation (7). 

As usual the inelastic processes includes both single-ladder exchange, as well as multiple 
interactions with a larger density of secondary partons. Analogous to the rescattering of a 
fast hadron in a heavy nucleus, we assume that the probability, w^, of events with parton 
multiplicity N times larger than that in a single ladder, is given by 



Unfortunately we cannot use this probability wn literally to describe the multiplicity distribu- 
tions of secondary hadrons. 

First, the distribution is affected by coherence effects. For example, in quark-quark scatter- 
ing, mediated by the exchange of t-channel gluons, the colour flow between the quarks cannot 
exceed the flow corresponding to octet (i.e. one-gluon) exchange. Thus, for low p^, when all A^ 
gluons act coherently, we expect the multiplicity of hadrons produced by hadronisation of the 
(A^-gluon) colour tube to be the same as that for one-gluon exchange, see also [37]. Secondly, a 
non-negligible fraction of the final hadrons may be produced via the fragmentation of minijets. 
These processes are beyond the 'pure soft' approach used in the present paper. Therefore below 
we concentrate on the processes with rapidity gaps. 



Wn = 




(31) 
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(a) SD 



(b) DD 



(c) DPE 




Figure 13: The cross sections of these three high-mass diffractive processes are calculated in this 
paper. The processes are (a) the dissociation of a single proton, (b) double dissociation (DD) 
and (c) central production (DPE). The shaded rectangles represent the interactions within the 
rapidity gaps, and are the result of complicated multi-Pomeron effects, see Fig. 11. 
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^ g-(f2fc+Oi)/2 
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Figure 14: The symbolic diagram for eq.(33) which describes single proton diffractive dissocia- 
tion process pp —>■ X + p. The origin of the various factors in (33) are indicated. 



4.4 High-mass diffractive dissociation 

We evaluate the cross sections for the three processes shown in Fig. 13, that is for the disso- 
ciation of a single proton (SD), for the dissociation of both incoming protons (DD) and for 
the central production of a system of mass M separated from the outgoing forward protons by 
rapidity gaps (DPE). In naive simplified models the latter process is called double- Pomeron- 
exchange. In reality the rapidity gaps shown as shaded 'rectangles' in Fig. 13 are the result of 
complicated multi-Pomeron effects, see Fig. 11. The results for the cross sections calculated in 
this paper correspond to rapidity gaps with Ay > 3. 



The cross section for single diffractive dissociation, ctsd 

We consider the diffractive dissociation of the beam particle i into a system of large mass M 
which occupies the rapidity interval from Y to y. It arises from the elastic scattering of an 
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/\ g-(nfe(yi)+f2i(y-yi))/2 



(l_e-§^^/2)2 



Figure 15: The symbolic diagram for eq.(36) which describes the double diffractive dissociation 
process pp —>■ X + X'. The origin of the various factors in (36) are indicated. 



intermediate parton c on the target /c, see Fig. 14. This elastic cross section, 

a = (1-6-^/2)2 (32) 

arises from the absorption of parton c as a consequence of the solution of the unitarity 
equation (2) for parton c, see (7). In this way we obtain a process with a rapidity gap in the 
interval {0,y). The resulting cross section for single diffractive dissociation is 

dy' 

with y' = Y — y = InM^. It is proportional to (i) the probability to find parton c in the 
interval dy', that is A e"'^^'' + ^'^^^ of (30); (ii) to the amplitude /«(?/') of the parton c - beam 
i interaction; (iii) to the gap survival factor S"^, that is the probability S'^ = e~-^^^^ to have no 
additional i — k rescattering; and (iv) to the elastic c — k cross section (1 — e~^/^)^. The opacity 
Q which drives this elastic amplitude is the opacity of the target, Qk{y)- 

The cross section for double diffractive dissociation, ctdd 

Double diffractive dissociation (DD) can be calculated in an analogous way. Here we have 
to consider the 'elastic' scattering of two intermediate partons, say parton c aX y = yi with 
y' = Y — yi and parton d at y = y2 < yi with y' = Y — y2. The corresponding opacity can be 
calculated as the difference between the beam opacities aX y = y2 and y = yi 

Q.(y - fe) - n.(Y - m) n,{Y ~ y^) 

The denominator, Qi{Y — yi), is the probability that an additional opacity 6Q is generated by 
the branch which contains the parton c. Alternatively the same quantity SQ can be calculated 
from the target side 

dil = - 1. (35) 
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Figure 16: Pomeron loop diagrams occurring in double diffractive dissociation. 



After the system of equations (29) and (30) is solved, both expressions give the same result, as 
may be checked numerically. 

The differential cross section for double dissociation is 



dyidy2 



^^_g-5n/2)2^g-(f7,(y2)+n,(y-s/2))/2^g-(c,(,i)m(y-j/i))/2 f.{Y-y^) fk{y2)S\ (36) 



The amplitude (1 — e"*^^/^) describes the exchange of Pomerons between partons c and d. In 
this way we allow for the Pomeron loop diagrams like those shown in Fig. 16. The factor 
l/n\ which arises in the expansion of the exponential e~^^/'^ accounts for the identity of the 
Pomerons arising in elastic c — d scattering. 

Double diffractive events, with a rapidity gap in the interval {y2,yi), may also come from 
two simultaneous single dissociation interactions: one with a gap in the interval {0,yi) and the 
other with a gap in the interval (?/2, Y)- This contribution to the DD cross section is given by 

,2 (SD*SD) 

= hiVi) 4(1/2) ^^ (37) 

dyidy2 

where 

/.(y) = /,(y') (l-e^^^J^)/")" Ae-(^^+^>)/2 (3g) 

Ikiy) = fk{y) (1 - e^^(^')/2)^ (39) 
are the probabilities of single dissociation (before the absorption), see (33); and y' = Y — y. 



The cross section for central DPE production, (Jdpe 

We can also calculate the cross section for the central DPE production of the process shown in 
Fig. 13(c). The central system has a mass given by = ^2 s, where = (1 — Xil) is the 
energy loss of the incoming protons. In comparison with the cross section of double dissociation 
of (36), where we consider the elastic scattering of the pair of intermediate partons c and d, 
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with rapidities yi and 1/2, now we have the inelastic c — d interaction accompanied by the elastic 
scattering of parton d on the 'target' k and parton c on the beam i. Thus the cross section is 
given by 

— — = t,k t,i b — — -, (4Uj 

"2/1 "2/2 Jk[y2) 

with 

Ej^ = (1 - e~^k{y2)/2Y ^ g-(nfc(j/2)+ni(y-2/2))/2 ^^-^-^ 
Ei = {1- e-^^«C^-yi)/2)2 ^ g-(Qj^(j;i)+Q,(y-yi))/2^ ^_^2) 

where the Ei and are the probabilities that the partons c and d participate in elastic 
interactions (multiplied by the probabilities to find these partons in the rapidity intervals dyi 
and dy2 respectively). The rapidities are Y — yi = — ln,^i and y2 = — In ,^2- The last factor in 
(40) describes the probability that the c — d interaction produces the central system of mass 
M; in other words the factor (1 — e^''^^^^/^)^ in (41), which corresponds to elastic scattering, 
replaces the factor fkiyi) in the whole amplitude. 



Asymptotic behaviour of diffract ive dissociation 

At very high energies, Y — \n s ^ 1, the whole irreducible amphtude given by (29) and (30) 
increases slowly 

/ ~ i In(yA), (43) 
see (27). On the other hand, the cross sections for single and double diffractive dissociation, 

asD ~ \ (AF)-i/\ (TDD ~ (AF)-i/\ (44) 

decrease as (Ay)"^/'^, due to the behaviour of the gap survival factor S'^ — ~ {AY)~^^^. 
In the product Ae~'^^'"'"^'=^/^ fi ~ l/Al^, the factor Y in the denominator is compensated by 
the size of the available rapidity interval, / dyi ~ Y. 



5 Accounting for the size of the coUiding hadrons 

Coming to the real 4-dimensional world, we may exploit the fact that the slope of the bare 
Pomeron trajectory needed to describe the elastic pp data {a' < 0.1 GeV~^) is very small in 
comparison with the slope of the proton form factor, Bq ~ 2.5 GeV~^. That is, the size of 
the Pomeron is much less than the size of the proton. Thus we may use the " heavy Pomeron 
approximation" [38], solving the system of equations (29) and (30) written in zero transverse 
dimensions, but with initial conditions Qkipk) and gi{bi), which depend on the position of an 
effective heavy Pomeron with respect to the centre of the beam and target protons in impact 
parameter space. 
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The pp impact parameter is b — b^ — bi, where here, and below, the vectors he in the 
2-dimensional transverse plane. In 6-space the input conditions for /j are given by the usual 
Fourier transform of the proton- Pomeron vertex /3j 

9i{bi) = J j^:^ Pi{q^) exp{iq-bi), (45) 

and a corresponding identical equation for g^. The whole irreducible amplitude fn-, which 
should be used to calculate the multichannel eikonal flik of Section 3.2, is given by 

fik{Y, b) ^ J d%d\ U{Y; 6,, h) S{h -%-b) (46) 

with 

fikiY, bi, 4) = fk{y = Y, bk, bi) Qiik) = fi{y' = Y", 6^, h) gkih)- (47) 

— * — * — # 

We keep both the arguments b^ and bi since, at each impact parameter b, we need to know the 
parton densities in both the beam and target protons. Indeed, the parton density fk{y) obtained 
from the evolution (29) of the target partons depends not only on the impact parameter inside 
the target proton, b^, but also that, 6j, for the beam proton through the external 'background' 
field fli{y') in (29). In terms of fik{Y, b), the total cross section takes the form 



= 2 1^^6(1 -e-^**(W)^ (48) 



see (3) and (6). This completes our discussion of the elastic amplitude fik- 

To calculate the high-mass contribution to o"sd we have just to replace fik in (46) by ex- 
pression (33) for the 'local' (that is fixed b) contribution to ctsd- The cross section takes the 
form 

^ ^ / f.^y') (1 _ e-^^/2)2 A e-("^ + S\b) d%d%i (49) 
dy' J 

with y' — Y — y — InM^. In terms of the (experimentally measured) longitudinal momentum 
fraction xl of the leading proton, we have dy' = dx^/il —xl). In (49), the eikonal gap survival 
factor S"^ depends on the pp impact factor b — bk — bi, and the dependence of /j, fk on the 6's 
is implicit. Note that now the factor S'^ accounts for the possibility of rescattering at any point 
in the 6-plane. That is 

S'{b) = exp(-4(r,6)). (50) 

The cross section of double diffractive dissociation is given by an analogous expression to 
(49), but in which now the integrals have integrands given by (36) and (37). However, in the 
c"dd *^°'* contribution of (37) we have to account for the possibility that the two simultaneous 
'single dissociations' li and Ik may take place at different impact parameters. That is we add 
the contribution 



d 



2 (SD*SD) 



a 



DD 



dyidy2 



j d%S\b)ii{y,,b)ik{y2M, (51) 
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where, in analogy to (46) , 



ij{y,b) ^ Jd%d%Ij{y-Xh)S{h-bi-b) (52) 

with j = i, k. The DPE cross section for central production is of the same form as (49) but 
with (40) as the integrand. 

We emphasize that everywhere in this section we are not dealing with the original proton 
beam and target, but rather with the Good- Walker diffractive eigenstates i,k = (j)i,4>k- Note 
that in this way we account for any low-mass intermediate states, that is we include all processes 
of the type pp — > X+p, pp ^ X + N*, .... Finally we take the sum over the different eigenstates 
as described in Section 3. 

In terms of Reggeon Field Theory, the exponential factors exp(— Q/2) create a very compli- 
cated system of multi-Pomeron exchange diagrams. However, as was mentioned in Section 3.1, 
actually the physical interpretation of such factors is simply the probability that no inelastic 
scattering occurs. Therefore, in terms of the parton approach, the equations generate only 
simple ladder diagrams, like Fig. 9 or Fig. 10, modified so that the probability (or splitting 
function) to produce a new parton (like c) is reduced by the 'survival probability' of the parton 
in the background fields. 

Now let us calculate the double differential diffractive cross section, d'^aso/dy'dq^, for — > 
X +p, where qr is the transverse momentum of the recoil proton. To do this we have to perform 
the Fourier transform of the pp ^ X + p amplitude with respect to the impact parameter b^, 
that is 



dy'dqj' 



j j \Jfi{y') (1 - e-^k/^Y ^ g-(n,+f2,)/2 ^2(g) ^iqT-hu ^2^^ 



47r ^ ' 



6 Practical details of the model 



Before we present the numerical results, it is convenient to list further details of the model. 

• The low-mass excitations are included in the Good- Walker formahsm (see Section 3). 
Therefore, to avoid double counting they should be omitted from the irreducible amplitude 
/ife. Hence, we introduce a threshold ?/o, and start (finish) the rapidity evolution of 
equations (29) and (30) at y = yo {y = Y — yo). We choose y^ = 2.3, a value which is 
often used in Regge Field Theoretic calculations, and which effectively accounts for the 
next-to-leading log corrections to the BFKL equation, see [39, 29]. Such a cut means that 
the low-mass dissociations with M < 2.5 GcV arc accounted for within the Good- Walker 
eikonal formalism, while larger masses are described in terms of multi-Pomeron diagrams. 



22 



• We consider two versions of the proton wave function, or, rather, of the 0j-Pomeron 
couphng Pi{t): 

(A) m = A(0) V(t) (54) 

(B) m = A(o) v(tm/{m)) (55) 

where in the first version, each component i has the same form factor, whereas in the 
second version the form factor of a component with a large cross section (large Pi{0)) 
is steeper than that for a component with a smaller cross section. In case (A) each 
component has the same transverse size and a larger cross section arises from a larger 
parton density. In version (B) the maximum parton density (which occurs at b—0) is the 
same for all components, and a larger coupling is caused by a larger transverse size. 

• The form factor of the 0j-Pomeron vertex is taken to have the form 

at 

(1 - t/aiY 

We take a ~ 0.1 GeV~^ in the exponent in the numerator. It is introduced solely to 
provide better numerical convergence of the Fourier transform (45). 

• The dispersion 7^ of the coupling, (15), is fixed by the experimental data on low-mass 
diffractive excitation. The analysis of the existing fixed target data require 4°D*^/^ei ^ 
0.3 [3, 40]. This value is consistent with the CERN ISR measurements of the excitations 
into particular channels {Ntt, Ntttt, AK, ...) with M < 2.5 GeV [41], and corresponds to^^ 

4°^^ ~ 2 mb at V5 = 31 GeV. (57) 

On the other hand in Ref. [42] the cross section extracted from measurements in the 
region M^/s < 0.01 was o"s*^*^ ~ 4 mb for s ~ 1000 GeV^. However in this experiment 
the momentum resolution AM^ was comparable to the whole size of the interval. 
Ax I = 0.01. Thus the quoted value, c^^^'^'^ ~ 4 mb, includes some contribution from 
masses larger than 2.5 GeV and, more important, may also contain some admixture of 
elastic events. In order to study the possible effect of a larger cross section than (57), 
we repeat the entire analysis with CTgo^ — 3 mb. (A value just below 3 mb was also 
measured in [43] .) Thus we will consider four versions of the model, which we denote by 
B2, B3, A2 and A3, where the letter A or B corresponds to the choice of couphngs given 
by (54) or (55), and the number 2 or 3 refers to crgo^ in mb for s ~ lOOOGeV^. 

• For both the n — 2 and n — 3 channel eikonals we take a simplified and extreme decom- 
position of the proton wave function 

\P) = -^t\<Pi)- (58) 

i=i 
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Here, and in what follows, the value of asu accounts for dissociation of both colliding particles. 
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• The value of A was defined in Section 4.2 {g'^ oc which leads to Q = A/). It controls 
the probability of high-mass diffractivc dissociation. A was tuned to be consistent with 
the measured CDF data [27, 28] for daso/dxLdt. 

• Prom the theoretical point of view the "heavy Pomeron approximation" is not totally 
satisfactory. To satisfy i-channel unitarity we need a non-zero slope of the Pomeron tra- 
jectory. At least, the two-pion loop contribution to Pomeron exchange should be included 
to satisfy unitarity at the nearest t = 4m^ branch point. However this contribution turns 
out to be small. We therefore set the slope a' = 0, which provides reasonably agreement 
with the data in the region of interest. Since we take a' = 0, we expect to underestimate 
the clastic slope, B^i at the LHC by up to 2%. Finally, the intercept of the bare Pomeron, 
q;jp(0), is left as a free parameter. 

7 Results of the model 

We fit the parameter oi of (56), which specifies the form factors of the Good- Walker eigenstates 
(f)i, and tune the value of a]p{0) = 1 -|- A, in order to describe the energy dependence of the 
measured total and elastic cross sections, and dad/dt in the CERN-ISR to Tevatron energy 
range. Recall that the parameters 7^ and A are specified by the available single diffractive 
data. A sample of the results are presented in Fig. 17^^ and Fig. 18. 
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A 
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0.53 


0.22 




0.9 


85 


(A2) 


0.40 


0.30 




0.42 


47 


(B3) 


0.65 


0.30 


1.80 


0.48 


38 


(B2) 


0.55 


0.33 


1.55 


0.275 


33 



Table 1: The values of the parameters in the various models; where (Jq is the square of the 
average value of the couphngs /3i(0): that is ctq = (/3i(0))^. By model (A2), for example, we 
mean that we use version (A) of the 0j-Pomeron couphng, (54), and we take mb, see 

(57). The reason for the missing entries for oi for model (A) is explained in the text. 

Even without any additional tuning of the parameters, the difference between the results of 
the two- and three-channel eikonal models is rather small (less than about 5 — 10%), once the 
dispersion 7^ of the couplings Pi{0) is fixed. So, below, we present results just for the three- 
channel case. The values of the parameters are listed in Table 1. Note that the predictions at 
the LHC energy are practically independent of the specific values of the vertices Pi{0) and of 

^^We thank Asher Gotsman for pointing out an incorrect labelling of the SD and DD curves in our previous 
version of the paper. 
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Vs (GeV) 

Figure 17: The description of the total and diffractive dissociation pp cross sections, obtained 
with a three-channel eikonal, for model (B2), shown by the continuous curves, and model (A3) 
shown by the dashed curves. The references to the data are as in [7] . Also shown by the curve 
denoted DL is the naive expectation for crtot obtained with a simple effective Pomeron pole 
(and secondary Regge contributions) [1]. The curve marked "fifajp only" is to demonstrate the 
importance of higher multi-Pomeron contributions, see the text. 
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do ,/dt (mb/GeV^) 




-t (GeV') 



Figure 18: The description of the elastic pp differential cross section for models (B2) and (A3) 
shown by continuous and dot-dashed curves respectively. The references to the data are as in 
[7]. 
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the chosen form of the t dependence of Pi{t), once we fix the value of low-mass dissociation, 
CTg^-^ (that is the dispersion 7^), and we satisfactorily describe the present data on the elastic 
cross section dagi/dt. Actually the t-behaviour of the vertex V(t) is strongly constrained by 
the data. If we choose a form different to (56), then after tuning the parameters to describe 
the data, the behaviour of V{t) turns out to be similar to that coming from (56). Also note 
that since the cross section for low-mass dissociation is screened by large rescattcring effects, 
the value of 7^ needed to describe the data is larger than the measured ratio C7gg"^^/2crei, see 
(16). In particular, we require a larger 7^ in models (A) than in models (B) to reproduce the 
same ratio 4°^^/2(7ei. 

The next observation is that it proved difficult to describe the duei/dt data using model (A), 
which has diffractive eigenstates of the same size. In fact we did not succeed with ansatz (56). 
More complicated forms were required. For model (A3) we took 

^^'^ ^ (1 - V0.58) ^ml-t ' ^'^^ 

while for model (A2) we have 

g0.05t 

where t is in units of GeV^. Even then, though the descriptions are reasonable, they are still 
not as good as those of model (B), see Fig. 18. 

In summary, the description of the total cross section in Fig. 17 is satisfactory, except at 
the lowest energies considered, where we have to add the secondary Regge contributions, which 
are outside our analysis. The differential elastic cross sections, shown in Fig. 18, are also 
well described. At the bottom of the plot we show the cross sections for the elastic scattering 
of two individual eigenstates 0i and 03 at the LHC energy. For 0i, which has a large cross 
section and large transverse size, wc sec that the diffractive dip occurs at a rather small t, 
—t ~ 0.2 GeV^. For the component 03, with the smallest cross section, there is no dip for 
— t < 1 GeV^. However, after we take the sum over all the (j)i(j)k contributions, we observe that 
the resulting amplitude is structureless. Note that the interference is important. After the 
dip the 0101 amplitude changes sign. Indeed, in model (B2), for —t > 0.5 GeV^ the whole pp 
daei/dt becomes smaller than that due to the small size component 03 alone. 

The best way to check the transverse structure of the proton is to consider diffractive deep 
inelastic scattering (DDIS), in particular the ^-dependence of elastic J ftp electroproduction on 
the proton. Here the beam particle corresponds to a 'heavy photon', which acts as a small 
size probe of the target. In our approach, we have simply to use a beam particle cigenstate 
with a Good- Walker wave function of small size and small absorptive cross section. We then 
obtain the t distribution shown in Fig. 19. If is sufficiently large, the transverse size of 
the incoming beam eigenfunction is much less than that of the target, so the t distribution is 
independent of Q^. Even though in model (B) we have components of quite different transverse 
size, for example in model (B2) 

Rl : Rl : Rl = 1.37 : 1 : 0.26, (61) 
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the resulting t distribution is very smooth. It is of the form 

~ e^* with 5 = 4-5 GeV-^ (62) 



dt 

in agreement with HI data [44]. 




Figure 19: The t distribution predicted for diffractive J/%Ij electroproduction using models (B2) 
and (A3). 



In Fig. 20(a) we present the cross section for high- mass single diffractive dissociation in 
the kinematical region measured by CDF [27, 28]. To describe the data at relatively large 
^ = {1 — xl) we add an Xi-independent secondary Regge contribution, RRP 

aa ^ aa^ ^ mb/GeV^. (63) 
dxidt axLttt 

The last term was normalized to the large ^ data and is almost negligible for ^ < 0.01. The 
energy dependence for fixed ^ = 0.01 is shown in Fig. 20(b). In Fig. 20(c) we show the 
t-dependence of this cross section. Again, despite the fact that model (B) has components of 
quite different radius, the t-dependence is well described by a single exponential with slope 
i? ~ 8 GeV~^. This is compatible with the CDF measurement [28]. The major difference 
between the predictions of models (A) and (B) is expected in the central diffractive cross section. 
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Figure 20: The behaviour of the cross sections for diffractive (a,b,c) single dissociation and (d) 
central production obtained for two of the models, namely (A3) and (B2). The four plots are 
(a) the description of the CDF data [27, 28] for single dissociation; (b) the energy dependence 
of dasn/dtdy, where y = — In^ with = {1 — xl); (c) the t-dependence of dasD/dtdx^, and 
(d) the y2 dependence of da-D^-^/ dyidy2 for = 0.05 and 0.005, corresponding, respectively, to 
proton taggers at 220 m and 420 m from the interaction point in the LHC experiments. 
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da^pE- This is illustrated in Fig. 20(d) for two values of the forward proton momentum fraction 
Xi = 1 — ^, which correspond to the acceptance of roman pots placed at 420 m = 0.005) and 
at 220 m (^ = 0.05) from the interaction point in an LHC experiment [45, 46]. The predicted 
values of the cross section have been obtained by integrating over the transverse momenta of 
both forward protons. 

The impact parameter shape of the elastic amplitude, Td, is shown in Fig. 21 for models 
(B2) and (B3). Since both models are tuned to describe the dad/dt data of Fig. 18 it is not 
surprising that the resulting amplitudes are close to each other. In fact, lmTci(6) should be 
regarded as an experimentally measured quantity, up to a small correction caused by the real 
part of the amplitude. Indeed, we have [47] 




I'-M^) = / 1/^ M^t-) (64) 



where = \t\ and p = ReTei/lmTei. From this point of view the latest GLM model is completely 
inconsistent with the data, as shown by the elastic amplitudes of Fig. 10 of the first reference 
in [9]. 

Note that the energy dependence of Tei is very weak in going from Tevatron to LHC energies. 
This is a common feature of models which have a larger number of diffractive eigenchannels and 
which include high-mass diffractive dissociation. First, if we have many different eigenstates, 
the part of the cross section which comes from the eigenstate with the smallest cross section 
reaches saturation much more slowly than the other components'^. Next, once we include the 
triple- (and multi-) Pomeron vertices we change the approach to saturation from the eikonal 
form, 

T = (65) 
to the power-like behaviour of the type 

see, for example, the Schwimmer formula given in (20). Here we have Q = = 

As we described in Section 3.3, the rapidity gap survival factor caused by multichannel 
eikonal rescattering depends on the structure of the matrix element Mik- In Fig. 22 we present 
the impact behaviour of the survival factor averaged over the diffractive eigenstates for two 
different assumptions about Aiik- First we assume that Aiik does not depend on the particular 
eigenstate, that is Aiik = M.- An example is the exclusive production of a Higgs boson via 77 
fusion. The probability to emit a photon is given just by the electric charge of the incoming 
state, and therefore is the same for each 0j. The corresponding curves are denoted by 77 — s>H 
on Fig. 22. The predicted 5*^(6) is approximately model independent. Another possibility is to 
assume that A4ik oc PiPk, which corresponds to Higgs production by Pomeron- Pomeron fusion. 
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This observation was also noted in Refs. [9] . 
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b (fm) 

Figure 21: The impact parameter shape of the elastic amphtude obtained from model (B) at 
the Tevatron and LHC energies. Two values, 2 and 3 mb, of o"g°j^^^ at s ~ 1000 GeV^, have 
been used to fix the parameters of the model. 



31 




0.5 1 1.5 2 2.5 3 3.5 



b (fm) 

Figure 22: The impact parameter behaviour of the rapidity gap survival factor S"^ averaged 
over the diffractive eigenstates 0j. Two assumptions are made about the matrix element for 
producing the system from eigenstates i, k, that is the amplitude J^ik in (18). First, ^Aik = M, 
for which Higgs production via 77 fusion is an example. Second, ^Aik oc PiPk, for which 
Pomeron-Pomeron fusion is an example. In each case the results are shown for four versions 
of the model: (A2), (A3), (B2) and (B3). Here the value shown for S'^{s,b) is equivalent to 
exp(— fi(s, 6)) in the case of a one-channel eikonal. It is the gap survival probability in an 
interaction with fixed proton-proton impact parameter b. To get the complete survival factor 
S'^, as given in either (17) or (18), we must also average S'^{s, b) over b with a weight given by the 
6-dependence of the matrix element Ai{b). For example, if, as in (69), A4{b) oc exp(— 6^/45) 
then the full ^ = ^ <fb W{s, b) exp(-6V2E)/27rE. 
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The corresponding curves are denoted IP IP — >H in this case. Here we have a smaller S'^{h) since 
the largest contribution comes from the diffractive component with the strongest absorption. 
Evidently the models with the larger o'gg'*^ give a smaller S'^(h) due to the larger contribution 
of the excited [N*) states to the absorptive correction. 

8 Conclusions 

We present a 'new generation' model for soft high energy hadron interactions, which includes a 
whole series of multi-Pomeron vertices. Of course the vertices corresponding to the interactions 
of many Pomerons are not known. Therefore such a model may contain an infinite number of 
parameters. To make progress we have assumed an eikonal-like form 

gl^ oc A"+"^-2 (67) 

for the m n Pomeron coupling. We emphasize that this assumption is much more reasonable 
than the assumption = for n + m > 3, that is keeping only the triple-Pomeron coupling 
g^ip. In fact, if wc were to neglect all the higher Pomeron vertices and keep only g^jp, then 
we would obtain much stronger absorption. Indeed, the curve denoted 'g^jp only' in Fig. 17 
was calculated with exactly the same values of the parameters as in the full model (B2), 
and demonstrates almost complete saturation starting already at the Tevatron energy. This 
illustrates the important role of the higher multi-Pomeron vertices. 

Another argument in favour of our approach is that it has a 'transparent' partonic interpre- 
tation. All the multi-Pomeron interactions are collected in the exponential factors e""*^^'"'"^*^'^^ 
which describe the absorption of intermediate partons during the evolution of the parton cas- 
cade, see (29) and (30). 

The main goal of this paper was to formulate a new and consistent approach to high energy 
soft hadronic interactions. However, it is informative to check that the approach is compatible 
with the present data. We found that already with a minimum number of parameters (see 
Table 1), all versions of the model can be tuned to give a reasonable description of the available 
data on soft hadronic interactions in the CERN ISR- Tevatron energy range, and to give similar 
predictions for the LHC energy. The reason is that the energy behaviour of the amplitudes is 
very smooth; there are no thresholds and everything depends on Ins, which does not change 
much in going from the Tevatron to the LHC energy. Nevertheless, an important prediction of 
the model is the relatively low value of the total cross section expected at the LHC: Utot — 90 
mb. This is an inescapable consequence of the absorptive corrections caused both by the low- 
mass intermediate states in eikonal rescattcring and the high-mass excitations described by 
multi-Pomeron interactions^^. The prediction for (Xtot is still compatible with the cosmic ray 
values [49], although it lies below the central cosmic ray expectations. 

^'^A similar value of o"tot was foimd in a simple two-parameter phenomenological analysis [48]. However high- 
mass diffraction is not embodied in the parametric form and much less diffractive dissociation was predicted. 
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Figure 23: The cross section dasB/dx^ for single dissociation integrated over t at the LHC 
energy resulting from four models: the continuous (dotted) curves are due to the B2 (B3) 
models, while the upper (red) dotted and dot-dashed curves are for models (A3) and (A2) 
respectively. For comparison we also show by a dashed (blue) curve the cross section obtained 
from model (B2) at the Tevatron. The secondary Regge contribution is included in the same 
way as in Fig. 20(a); it is relatively very small for (1 — xl) < 10~^. 
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Another example of very flat energy behaviour is that shown for daso/dxL in Fig. 23. We 
see that this single dissociation cross section changes very little^^ in going from the Tevatron 
(dashed curve) to the LHC (solid curve for model (B2)). This quantity, which is interesting 
in its own right, is of practical importance since it governs the rate of the so-called pile- up 
background to central exclusive production processes which are caused by the overlap of two 
soft diffractive events simultaneously with a hard scale event [23, 50]. 

Note also that the rapidity gap survival factor also decreases slowly with energy. For exam- 
ple, we can see from the final entries of Table 2 that 

(Tevatron) / ^^^(LHC) -1.5. (68) 

This means that the present experiments at the Tevatron can be used to check the theoretical 
expectations of the models at the LHC energy. 

Actually the survival factor, which governs the rate of events in a particular diffractive 
process, must be calculated by averaging the value of S'^ip) with the forms of the matrix 
elements M.ik{h), as in (18). The values presented in Table 2 assume, for simphcity, that^^ 

Mik oc exp(-6V4S). (69) 

(Of course, in general, each amplitude M.ik may have its own slope -Bjfc.) As examples, we 
choose the values B — A and 5.5 GeV~^ to be the same as in [7]. After we tune the model to 
the soft data, we find that we still obtain almost the same survival factors at the LHC energy; 
compare the entries in Table 2 with those in the "CD(FPS)" column of Table 1 of [7]. 

A sample of the results for the total and the diffractive dissociation cross sections were 
presented in Fig. 17, and for dae\/dt in Fig. 18. In Table 2 we also list the values of the cross 
sections for models (B2) and (B3) at three energies, namely ^/s = 1.8 GeV, 14 TeV and 10^ 
GeV. Note that model (B) is favoured, for reasons explained below. We see that with increasing 
energy the cross sections violate the Pumplin bound [51] 

O-el + CTSD + CTdB < CTtotf^- (70) 

However this bound is only justified for low-mass dissociation described in terms of Cood- 
Walker diffractive eigenstates. (Actually the bound was originally proved for dissociation on 
nuclei.) For the low-mass components of the diffractive dissociation cross sections we see that 
there is no contradiction with the bound. 

An important ingredient of the analyses of soft high energy hadron interactions are the 
data on diffractive dissociation. The data that are available at present are fragmentary. No 
experiment has covered the whole kinematic range of t and rapidity. Moreover, at low M and 

^^The growth of the cross section caused by the Pomeron interactions (A > 0) is compensated by stronger 
absorption, that is by the decrease of the gap survival factor 5^. 

^■'Notc that Ref. [7] denotes the slope as 26 = B, where here b is not to be confused with the impact 
parameter. 
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Tevatron 


LHC 


= 10^ GeV 


O'tot 
0"el 


74.0 (73.9) 
16.3 (15.1) 


88.0 (86.3) 

20.1 (18.1) 


98.0 (94.3) 
22.9 (20.0) 


C^SD 

_lowM 
^SD 
hiffhM 


10.9 (12.7) 
4.3 (6.0) 
6.5 (6.7) 


13.3 (16.1) 
5.1 (7.0) 
8.1 (9.1) 


15.7 (17.7) 
5.7 (7.9) 
10.0 (9.8) 


_lowM 
highM 

(highM*lowAf) 
(SD*SD) 


7.2 (8.7) 
0.2 (0.5) 
4.5 (4.0) 
2.1 (3.6) 
0.4 (0.7) 


13.4 (12.9) 
0.2 (0.5) 

9.3 (5.9) 
2.9 (5.2) 
1.0 (1.3) 


17.3 (21.1) 
0.2 (0.6) 

11.7 (12.9) 
3.8 (6.0) 
1.6 (1.6) 


(5 = 4) 
(5 = 5.5) 


0.027 (0.018) 
0.048 (0.032) 


0.017 (0.012) 
0.032 (0.023) 


0.013 (0.009) 
0.025 (0.018) 



Table 2: The cross sections (in mb) predicted at the three energies using model (B2). Also 
shown in brackets are the predictions for model (B3). The various components of the single 
and double dissociation cross sections are also listed. The final entries are the survival factors 
of the rapidity gaps for exclusive diffractive central production, pp ^ p + X + p, for two values 
of the slope 5 = 4 and 5.5 GeV"^. The expression for a§^*^^^ is given in (37). 

low t it is hard to avoid contamination from elastic events. For these reasons we prefer to 
compare the theoretical predictions with measurements of differential cross sections, rfcTg^^/ dt 
and (Pag^^^ /dtdM'^ , in regions where they were actually observed, rather than with integrated 
cross sections, ctsd and ctdd, which were obtained by extrapolating the data into unmeasured 
regions using some simplified model. 

In order to further constrain the parameters of the model so as to obtain more precise 
predictions, it is important to make accurate measurements at collider energies, of the single 
diffractive dissociation cross section for low masses a^^^^^ and of its t dependence da^^^^ /dt, 
and of central diffractive production daj:,pE/dyidy2. The latter two measurements will provide 
information on the transverse sizes of the various diffractive eigenstates </){. Such measurements 
can be performed at the LHC^^, for example, in the TOTEM experiment [52]. 

We have considered two different models for the 0j-Pomeron couplings: model (A) in which 
all the components have the same size, and model (B) in which the sizes differ, see (54) and 
(55). There are three observations which appear to favour model (B). First, model (B) gives a 
somewhat better detailed description of the various data for daei/dt. Indeed, recall for model 
(A3) we were required to choose a complicated form of the residue V{t) to avoid a contra- 
diction with the data. Second, Model (B) predicts itdd which is better agreement with CDF 

^*We thank Risto Orava for discussion of this issue. 
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(a) (b) 

Figure 24: Typical final states which are included in the model: (a) diagrams from multi- 
channel eikonal rescattering, (b) diagrams with and without a rapidity gap from enhanced 
multi-Pomeron interactions. 

measurements^^. Third, model (A) requires a very large dispersion 7^ of the Pomeron-0j cou- 
plings, which tends to be incompatible with the data for soft reactions at fixed target energies, 
see [40]. Nevertheless, at present we cannot completely reject model (A), and so we have also 
presented the results of this model to indicate the possible spread of the predictions. 

In all the models the value of A was about 1/3, in agreement with the evaluation of g^ip 
based on the HERA data for 7p J/ip X [10]. 

In our description of high energy 'soft' hadron-hadron interactions we need to start with a 
rather large 'bare Pomeron' intercept: a]p{Q) = 1 + A with A ~ 0.5. Note that this is con- 
siderably larger than the intercept A ~ 0.08 of the naive effective Donnachie-Landshoff-type 
of Pomeron [1]. The reason is that it is necessary to compensate the effects of the screening 
corrections originating from the enhanced multi-Pomeron diagrams. Recall that in the present 
model we use the heavy Pomeron approximation [38], that is we set the slope of the Pomeron 
trajectory a' = 0. Nevertheless this extreme choice is still consistent with the data. The 
shrinkage of the elastic differential cross section data of Fig. 18 is reproduced in our approach 
by a stronger absorptive effect at low impact parameters with increasing energy. Our prelim- 
inary estimates indicate that a small value of a'{a' < 0.05 GeV"^) may possibly improve the 
description of the data^^, but we find a larger a' {a' > 0.05 GeV~^) appears to be definitely 
ruled out, since the main shrinkage of the elastic peak comes from absorptive effects. 

In summary, we have been able, for the first time, to present a theoretically fully consistent 
description of high energy soft hadron interactions. The procedure incorporates both a multi- 
channel eikonal and multi-Pomeron interactions. In terms of Reggeon Field Theory, it sums up 
the eikonal diagrams in Fig. 2, together with the enhanced multi-Pomeron diagrams like those 

CDF [53] have measured the cross section for double dissociation with a rapidity gap Ay > 3 enclosing y = 
and find Actdd = 4.4 ± 1.2 mb at 1.8 TeV and 3.4 ± 1.1 mb at 630 GeV, as compared to the predictions of 
Actdd =4.4 mb and 3.7 mb, respectively, of model (B2). Here Actdd denotes the portion of the cross section 
for double dissociation coming from the kinematic region where it was actually measured. 
^^Simultaneously, the inclusion of a non-zero a' will result in a smaller value of A. 
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in Figs. 5 and 16. In terms of final states these diagrams generate processes with different 
densities of secondary particles. The eikonal diagrams generate processes like those shown in 
Fig. 24(a), while the enhanced multi-Pomeron diagrams may generate processes with rapidity 
gaps, as well as those with different densities of secondaries in particular rapidity intervals as 
shown in Fig. 24(b); see the discussion in Section 4.3. 
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